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WIDE-BAND  SIGNAL  LEVELS 
IN  SONAR  MODELING 

I.  INTRODUCTION 

(U)  Whenever  possible,  narrow-band  (spectral)  mathematical  models 

are  used  for  predicting  signal  levels  in  sonar  systems;  occasionally, 

however,  models  are  required  for  systems  whose  bandwidth  is  so  large 

that  a narrow-band  approximation  cannot  be  justified.  In  such  cases, 

the  spectral  levels  must  be  integrated  over  the  frequenc-  band  of 

interest.  The  integration  is  commonly  performed  using  numerical 
1 2 

techniques,  which  are  relatively  time-consuming  ant  inaccurate. 

An  analytic  expression,  which  completely  eliminates  the  need  for  using 
the  less  efficient  numerical  techniques,  has  been  derived  ami  is 
presented  in  section  II.  It  is  applicable  whenever  the  spectral 
intensity  can  be  approximated  as  a series  of  straight-line  segments 
when  plotted  as  a function  of  frequency  using  log-log  scales.  This 
is  very  often  the  case  in  practice,  so  that  the  results  presented 
here  have  wide  applicability  to  practical  sonar  models. 

I I . BAND-LEVEL  FOR  A 1 'IDE-BAND  SIGNAL 

(U)  The  total  intensity  for  a wide-band  signal  can  be  obtained  by 

integrating  the  spectral  intensity  i over  the  frequency  hand  of 
a 

interest  [f  , f ] 

f6 

i = / i (f)  df,  (1) 

B ' s 

a 
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where 


iD  is  in  units-of-power/cm2* 
b 


ig  is  in  units-of-power/cm2/Hz ; 


f is  in  Hz. 


(U)  When  referring  to  intensities  and  intensity  ratios,  for 
convenience  we  shall  use  lower-case  symbols  for  the  actual 
intensity  and  upper-case  symbols  to  indicate  the  same  quantity 
expressed  in  decibels  with  respect  to  1 ubar  (dBM)  or  decibels 


(U)  Assume  that  i can  be  approximated  over  the  band  .if  interest 
using  an  expression  of  the  general  form: 


is(f)  = bf“ 


Then  I (f)  in  dBM  can  be  written: 
s 


Ig(f)  = 10  m log ( f ) + B 


which  is  the  equation  of  a straight  line  if  log(f)  is  the 
variable.  This  form  occurs  frequently  in  sonar  systems 
analysis  when  data  points  are  plotted  in  dBM  on  semi-log  paper. 


*The  reference  power  level  is  approximately  0.64x10  watts/cm2 
for  water3. 


isar.i 
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Substituting  eqn.  (2a)  into  eqn.  (1)  and  integrating  yields: 


(m+1) 


,,m+l  fm+l. 

a ;; 


m 4 -1, 


(3a) 


or 


(U)  Equation  (3)  solves  the  problem  completely,  in  principle.  Its 
application  to  practical  situations  arising  in  sonar  modeling 
illustrated  in  sections  III  and  IV  are  thought  to  be  both  useful  and 
novel. 

(U)  It  should  be  pointed  out  that  the  intensity  expression  given 
in  eqn.  (2a)  can  be  broken  down  into  a product  of  n similar  factors, 
each  corresponding  to  an  individual  sub-system  of  the  total  sonar 
system,  and  eqn.  (3)  is  still  applicable.  In  this  case,  in  eqn.  (2b), 


and 


(4a) 


* ■ I'. 

i=l 


(4b) 
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III.  APPLICATION  TO  ORDERED-PAIR  DATA 

(U)  Frequently,  spectral  intensities  Ij  are  measured  or  calculated 

only  at  a discrete  set  of  frequencies  (f^,  f2,...,f  ) in  the  frequency 

band  of  interest.  Such  information,  when  grouped  in  data-pairs 

(frequency,  intensity)  and  arranged  in  increasing  order  of  frequency, 

is  called  ordered-pair  data:  (f,,  I,).  (f„,  1„) (f  , I ). 

-L  i L L n n 

(U)  When  intensities  are  available  only  at  a discrete  set  of 

frequencies,  the  approximate  intensity  at  any  included  frequency  is 

*1  2 

often  obtained  by  linear  interpolation.  * This  defines  a continuous 
functional  relationship,  consisting  of  straight-line  segments  between 
adjacent  data-pairs,  over  the  frequency  band  of  Interest,  called  a 
"piece-wise  linear"  function.  Under  these  conditions,  eqn.  f2b)  Is 
applicable  between  each  frequency-pair,  since  the  data  are  plotted 
using  a logarithmic  frequency  scale.  Of  course,  in  general  there 
will  be  n-1  distinct  m^'s  and  Bj’s.  Figure  1 will  clarify  this 
situation. 

(U)  It  is  clear  that  eqn.  (3)  applies  individually  to  each  linear- 
interpolation  region,  so  that  the  total  band  integral  for  ordered-pair 
data  can  be  obtained  as  a sum  of  n-1  integrals,  one  for  each  region: 

* "Linear"  interpolation  is  actually  "Logarithmic"  since  the  independent 
variable  is  assumed  to  be  log(f),  not  f itself. 


SCL&SSR' 


n-1  f . 


2 Jv, 

3-1  f. 


(f)df  , 


where 


Igj(f)  = 1-0  mj  log(f)  + B ; j = 1 n-1. 


Each  of  the  n-1  integrals  in  eqn.  (5)  is  evaluated  individually  using 
eqn.  (3a)  or  (3b),  as  appropriate.  The  values  of  m and  B are  not 

j j 

explicitly-  available  in  the  ordered-pair  data  given,  but  they  can  be 
calculated  using  expressions  derived  from  the  geometry  of  fig.  1: 


m = h+i1  j 

j 10  log/f1+l 


log  (f.) 


B = T — J — 

j j log/Vl 


(Ij+rIj) 


In  order  to  clarify  the  method  for  obtaining  the  total  band-level 
intensity,  a step-by-step  summary  follows:  (a)  Eqn.  (6)  and  the 
given  data  are  used  to  obtain  n-1  values  each  of  m^  and  . (b)  The 
band-level  integral  for  each  of  the  n-1  regions  is  individually 
evaluated,  using  the  values  of  m_.  and  for  that  particular  region, 
and  (3a)  or  (3b)  as  appropriate,  (c)  The  total  band-level  intensity 
is  obtained  as  a sum  of  the  individual  integrals  for  the  n-1  regions, 


as  indicated  in  equation  (5). 
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(U)  This  proposed  method  is  superior  to  the  numerical-integration 
procedures  commonly  used,  for  at  least  two  reasons:  First,  it  is 
more  nearly  mathematically  correct,  since  the  integral  in  each 
region  is  exact,  not  approximate.  (The  only  approximation  is  the 
linear-interpolation  assumption,  which  applies  to  both  methods.) 

Second,  it  is  virtually  always  much  more  efficient,  since  many  more 
regions  will  be  required  in  the  numerical-integration  procedure. 

(U)  In  addition  to  the  obvious  case  of  measured  signal-level  data, 
at  least  one  other  case  of  practical  importance  exists  in  which  the 
signal  levels  are  available  as  ordered-pair  data:  In  modeling  a wide- 
band passive  sonar  system,  the  signal  level  is  approximated  by 
predicting  the  effects  of  the  ocean  and  the  sonar  system  on  the  target 
noise  spectrum.  While  the  frequency  characteristics  of  the  sonar 
system  are  usually  fairly  well  known,  those  of  the  propagating  medium 
and  the  target  noise  spectrum  are  usually  available  only  in  a highly 
complicated  form,  or  as  measured  data.  In  either  case,  it  is 
reasonable  to  evaluate  the  spectral  signal  level  only  at  a discrete 
set  of  frequencies  in  the  band  of  interest,  In  which  case  the  method 
presented  here  for  obtaining  total  band-level  is  directly  applicable. 

IV.  APPLICATION  TO  INCOHERENTLY-SUMMED  LEVELS 
(U)  Frequently  two  or  more  signals,  each  of  which  can  be  expressed 
in  the  form  of  eqn.  (2),  are  incoherently  summed  to  obtain  the 
spectral  level  of  interest.  In  this  case,  the  very  simple  method  of 
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summing  the  m^'s  and  B^'s,  given  in  eqn.  (4),  is  not  applicable.  A 
general  expression  applicable  to  this  situation  is  derived  in  section 
IV. B;  a special  case  of  practical  importance  will  first  be  treated 
in  section  IV. A. 

A.  Empirical  Expressions  for  Sea  and  Ship  Noise 

(C)  Three  empirical  expressions  for  predicting  ambient  sea  and 

4 5 

ship  noise  are  available.  * The  equations,  which  are  based  on  the 
Knudsen  curves  and  CDR  Fridge's  data  are: 


Nq  = -42  -16.67  log  (F)  + 10.35  log 
Nw  = -42.9  -16.67  log  (F)  + 2h 
Nv  - -44  -16.67  log  (F)  + 1.5v 


(h) 

dBM/Hz 

(7) 

dBM/Hz 

(8) 

dBM/Hz 

(9) 

where  Nq,  N^,  and  Ny  are,  respectively,  the  ambient  ocean  noise,  the 
wave  component  of  ship  noise,  and  the  velocity  component  of  ship  noise; 

F is  the  frequency  in  kHz; 
h is  the  wave-height  in  ft; 
v is  the  ship  speed  in  kts. 

(U)  In  sonar  predictions,  the  noise  sources  whose  levels  are  predicted 
by  equations  (7),  (8),  and  (9)  are  assumed  to  add  incoherently.  Thus, 
eqn.  (7)  alone  predicts  the  ambient  ocean  noise  in  the  absence  of  the 
ship,  the  incoherent  sum  of  equations  (8)  and  (9)  predicts  the  total 
additional  noise  N^  due  to  the  presence  of  the  ship,  and  the  incoherent 
sum  of  all  three  equations  predicts  the  total  ambient  ocean  noise  level 


with  the  ship  present  N^. 
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In  equation  form, 

N = 10  log(n  + n ) dBM/Hz  (10) 

U W V 

and 

N = 10  log(n  + n + n ) dBM/Hz  (11) 

1 O W V 


(C)  The  expression  for  the  total  spectral  noise  level  obtained  by 
combining  equations  (7),  (8),  and  (9),  using  (11),  is  useful  in  itself; 
it  will  be  obtained  prior  to  applying  eqn.  (3)  for  the  total  band 
level.  Summing  the  intensities  incoherently,  we  have: 


nT 

or 


(10-4,,2-h1,035  + 10"4,29-10°,2h  + 10~A,4.10°,15v)  F-1’667 


(12a) 


N„ 


-40  + 10  log 


1.035 


10 


0. 2h 


10 


0. 15v 


1.5849  1.9498  2.5119 


-16.67  log (F)  dBM/Hz 


(12b) 


(C)  Since  eqn.  (12b)  is  in  the  proper  form,  eqn.  (3a)  can  be 
directly  applied  to  obtain  the  ambient  noise  band  level  in  the  water, 
for  any  frequency  band  [F^  , F£]  over  which  equations  (7),  (8) , and  (9) 
are  applicable: 


nb  = 


-10  + 10  log 


1.035 


10 


.0 . 2h 


10 


,0. 15v 


1.5849  1.9498  2.5119 


+ 10 


-0.667 


-0.667 


log 


0.667 


dBM 


(13) 
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where  30  d”  has  been  added  to  allow  for  integration  with  respect 
to  frequency  in  kHz.  The  argument  of  the  logarithm  in  the  third 
term  of  eqn.  (13)  will  be  positive  for  the  condition  F2  > F^ , as 
was  assumed.  The  denominator  in  this  term  can  be  combined  with 
the  other  constants;  it  has  a value  of  approximately  1.75874  dB. 

(U)  Equation  (13)  gives  the  noise  band  level  in  the  water.  A 
more  useful  level  is  at  the  output  of  a sonar  system.  Equations 
for  the  system  functions,  if  available  in  the  proper  form,  can  be 
combined  with  eqn.  (12b)  using  eqn.  (4).  This  has  been  dot  e for 
the  SQQ-23  PAIR  system,  and  the  results  will  be  given  here  as  &u 
example. 

(C)  The  PAIR  system  functions  required  to  obtain  the  output  noise 

level  are  the  system  directivity  index  and  amplifier  gain  G. 

2 

The  expressions  are: 


D = 11  + 20  log  (F) 

dB 

(14a) 

G =16.6  log (F) 

dB 

(14b) 

(C)  Equation  (32b)  and  (14)  can  be  combined  using  eqn.  (4).  It 
should  be  noted  that  quieting  terms,  amounting  to  1 dB  and  3 dB, 
have  been  applied  to  equations  (7)  and  (9),  respectively,  so  that 
the  constants  differ  slightly  from  those  previously  given.  The 
spectral  level  (integrand)  obtained  is: 


10 


CONFiOLm.rt; 


> i IL/L^  l 


b 


n (?)  = 10"1,1  • (10-4,3h1,035  + 10_4,29.10°,2h 
sp 


+ IQ"4’ 7 .10°‘15v)  F-2 


(15a) 


or 

h1-035  , 10°-2h  . 10°-1Sv’ 

1-9953  1.9498  5.0119 

- 20  log(F)  dBM/Hz  (15b) 


N (F)  = -51  + 10  log 
sp 


The  expression  for  the  corresponding  band  level,  over  a 
frequency  band  [F^,  is : 


= -21  + 10  log 


h1’035  io°*2h  io°-L5v 

1.9953  + 1.9498  + 5.0119 


+ 10  log 


dBM 


(16) 


where  30  dB  has  again  been  added  to  allow  integration  with  respect 
to  frequency  in  kHz.  In  the  FAIR  system,  F^  and  F ^ are  known 
constants,  so  that  the  last  term  in  eqn.  (16)  can  be  evaluated  and 
written  into  the  equation.  Actually,  there  are  three  different 
passive  bands  used  in  the  PAIR  system:  1.0  to  1.8  kHz;  1.0  to  2.5 
kHz;  and  1.0  to  6.0  kHz.  The  corresponding  constants  for  the  third 
term  in  eqn.  (16)  are:  -3.5218,  -2.2185,  and  -0.7918. 

B.  General  Expressions  for  Incoherently-Summed  Sources 
(U)  In  the  general  case,  there  will  be  n signal  sources  whose 
intensities  are  assumed  to  add  incoherently.  Suppose  that  the 
intensity  of  each  source  is  given  by  an  expression  of  the  form  of 
eqn.  (2b): 


11 
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Sk  = 10mic1°8(f)  + Bk‘»  k ■ !»•••!  n. 


The  corresponding  net  spectral  level,  assuming  incoherent-source 
summing,  is i 

k-1 


(U.)  Equation  (17)  can  be  integrated  using  eqn.  (3a)  or  (3b)  on 
each  term  as  appropriate.  The  resulting  intensity  over  a frequency- 
band  [f^  , f2]  is: 


(17) 


(18) 


V.  CONCLUSIONS 

(U)  Well-known  integral  relationships  have  been  applied  to  practical 
sonar  intensity  functions  to  obtain  equation  (3) , a general  expression 
for  total  wide-band  intensity.  The  basic  assumption  is  that  the 
spectral  levels  of  interest  can  be  expressed  as  in  eqn.  (2b),  which  is 
frequently  the  case  in  practice, 
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(U)  Equation  (3)  was  applied  to  the  important  special  cases  of 
measured  data  and  incoherently-summed  sources,  with  the  results 
given  in  sections  III  and  IV.  As  an  example  of  the  latter  case, 
expressions  were  derived  using  empirical  expressions  for  Knudsen's 
curves  and  CDR  Fridge's  data  for  ambient  sea  and  ship  noise. 
Results  for  the  passive  SQQ-23  PAIR  system'  were  also  presented. 

(U)  Computation  time  using  the  expressions  derived  will  normally 
be  very  much  less  than  that  of  previously-used  numerical  methods, 
yet  the  results  obtained  will  be  more  accurate. 
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